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We study condensation transitions in the steady state of a zero-range process 
with two species of particles. The steady state is exactly soluble — it is given 
by a factorised form provided the dynamics satisfy certain constraints — and we 
exploit this to derive the phase diagram for a quite general choice of dynamics. This 
phase diagram contains a variety of new mechanisms of condensate formation, and 
a novel phase in which the condensate of one of the particle species is sustained by 
a 'weak' condensate of particles of the other species. We also demonstrate how a 
single particle of one of the species (which plays the role of a defect particle) can 
induce Bose-Einstein condensation above a critical density of particles of the other 
species. 
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I. INTRODUCTION 



Condensation phenomena are observed in a variety of contexts. For instance, microscopic 
dynamics including particle diffusion, aggregation to form particle clusters, and fragmen- 
tation of these clusters, can be used to model a number of physical systems . In one 

n 

dimension, such models have been analysed within mean field theory |1] or using a scaling 
approach [2] to infer the existence of transitions between a fluid phase a condensate phase. 
Another context is the modelling of granular and traffic flow One such model is the 'Bus 
Route model' 0|, where there is a crossover between a regime in which the average velocity 
of buses is determined by the velocity of the slowest bus, and, above a critical density of 
buses, a regime in which the average bus velocity is limited instead by the high density of 
traffic. This crossover can be understood in terms of a condensation process. Further, several 
models have been shown to undergo phase separation in one dimension P, 0]. This phase 
separation can also be related to a condensation mechanism. In particular, a general crite- 
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rion has been proposed, predicting the existence of phase separation in 1-d driven systems, 
which appears to be widely apphcable This apphcabihty follows from the robust nature 
of the physical mechanism underlying the phase separation. This mechanism, which is also 
the generic mechanism for the aforementioned condensation phenomena, is understood in 
terms of condensation transitions in the zero-range process 0|. 

The zero-range process is a system of many interacting particles which move on a lattice 
— particles hop to adjacent lattice sites with hop rates determined by the number of particles 
present at the departure site. It provides insight into the behaviour of more complicated 
models because it is exactly soluble: the steady state assumes a simple, factorised form. 
Thus the condensation transitions, whereby a finite fraction of particles occupy a single 
site, are amenable to exact analysis. Condensation transitions in the single species zero- 
range process proceed through one of two mechanisms: (i) if the particle hop rates are site 
dependent, then above a critical density a condensate forms at the site where the hop rate is 
slowest — this mechanism is closely related to Bose condensation, (ii) if the particle hop rates 
depend on the number of particles present at the departure site, then above a critical density, 
and provided the asymptotic dependence decays to a constant value sufficiently quickly, a 
condensate forms at a site located at random — thus the transition is accompanied by a 
spontaneously broken symmetry. The former mechanism here demonstrates how disorder 
can induce condensation, analysis of the latter leads to an understanding of the generic 
mechanism of condensate formation applicable to the examples described above. 

A question which naturally arises then is what other generic mechanisms of condensate 
formation exist. To this end, we consider the generalisation of the zero-range process to 
two species of particles. We investigate how the interaction of the two species allows new 
mechanisms of condensate formation. 

The zero-range process with two species of particles was introduced in P, [l^ where 
the steady state was obtained exactly and shown to be given by a simple factorised form 

n 

provided the dynamics satisfy certain constraints. In [9], this was used to demonstrate a 



new mechanism of condensation transition for a specific choice of dynamics, and in the 
hydrodynamics were derived. Here, we show how the two species model may undergo a wide 
variety of condensation transitions, and we derive the phase diagram for the model for a quite 
general choice of dynamics. Also, we show that the steady state of the model can be mapped 
on to the steady state of the AHR model — a model which undergoes a transition between 
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a disordered (fluid) phase and a phase separated (condensate) phase. Thus, as in the single 
species model, the two species zero-range process exhibits transitions of a robust nature 
which can provide insight into condensation mechanisms in more complicated models. In 
particular, one perspective of the two species model is to consider one species as providing an 
evolving landscape upon which the other species in turn evolves. This evolution is coupled, 
thus condensation transitions are induced by the evol ving disordered background. Again, 
such interplay arises in a variety of physical settings . Our aim is to explore 

how this interplay can lead to novel condensation transitions in the two species zero-range 
process. 

We begin by reviewing in Section 2 the key equations of the steady state solution, which 
form the basis of the subsequent analysis. In Section IIIH we show how a defect particle 
(i.e. a single particle of one of the species) can induce a condensation transition in the 
particles of the other species. We consider in Section |TV| the case where the hop rates of one 
of the particle species depend only on the number of particles of the other species at the 
departure site. We show how to derive the phase diagram for this case and find that three 
distinct condensate phases can arise; numerical iterations of an exact recursion relation for 
the partition function yield results consistent with the predicted phases. In Section |3 we 
present a mapping between the steady state of the two species zero range process and the 
AHR model. We conclude in Section IVTl 



II. STEADY STATE 



We define the two species zero-range process on a lattice containing L sites and with 
periodic boundary conditions. On this lattice, there are particles of species A and M 
particles of species B. Particles of both species hop to the nearest neighbour site to the 
right, species A with rate u{ni,mi) and species B with rate v{ni,mi), where site / is the 
departure site and contains rii particles of species A and mi particles of species B . 

Since the steady state has already been derived in detail elsewhere 0, [l^ , we quote only 
the key results here. We define P{{ni}; {rrii}) to be the probability of finding the system 
in the configuration ({^/}; {rrii}), where {n/} = rii, . . . ,nL and {mi} = mi, . . . , mi. This is 
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given by a factorised form 

L 

P{{ni\, {mi}) = Zll^j^j JJ fijii, mi) , (1) 

1=1 

where Zl,n,m is a normalisation. The steady state ((T)) satisfies the steady state master 
equation if the factors /(n/, mi) satisfy 

u{ni,mi)f{ni,mi) ^ ^ v{ni,mi)f{ni,mi) ^ ^ 

f{ni-l,mi) f{ni,mi-l) 

The solution to these equations is 

f{ni,mi) = l[[uit,mi)]-'ll[v{0,j)]-' , (3) 

i=i i=i 



provided the hop rates satisfy the constraint 

u{ni,mi) v{ni,mi) 



(4) 



u{ni,mi-l) v{ni-l,mi) ' 

for ni,mi ^ — the choices of u{ni,0) and v{0,mi) remain unconstrained. We emphasise 
that instead of specifying the hop rates directly, we have the freedom to choose any desired 
form for f{n, m) and that we can then infer the hop rates from 

The normalisation Zl^n,m, defined in dT)), plays a role analogous to the canonical partition 
function of equilibrium statistical mechanics, and is given by 

L L L 

Zl,n,m= SiY.'^i-N)SiY,m-M)l[f{ni,mi) , (5) 

{ni},{mi} 1=1 1=1 1=1 

where the delta-functions ensure that the system contains the correct numbers of particles of 
each species. By writing the delta- functions in an integral representation, Zl^n,m becomes 

dz r dy [F{z,y)]^ 



'L,N,M 



2m J 27ri z^+iy^-^+i 
where the generating function, F{z, y), has been defined as 



(6) 



oo oo 



F(z,y) = ^ J]^W(n,m). (7) 

n=0 m=0 

We consider in Section HVl the limit L, N, M ^ oo, where pa = N/L and pb = M/L — the 
particle densities of species A and B respectively — are held fixed. In this limit, we assume 
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that the integral in equation © is dominated by the saddle point. The equations for the 
saddle point are 

d d 
Pa = z—lnF{z, y) , pb = y—lnF{z, y) . (8) 

Assuming the saddle point is valid, equations (jHI) determine pA and in terms of z and y 
and this amounts to working in a grand canonical ensemble. We note that for the saddle 
point to be valid, z and y cannot exceed the radii of convergence of F{z,y), since we must 
be able to perform the sum ((Zj) in the first place. Further, since all derivatives of F{z, y) are 
positive, the saddle point, if valid, must be unique. In Section HVl we will find that it is not 
always possible to solve the saddle point equations for all values of pa and in the allowed 
ranges of z and y. This phenomenon corresponds to a condensation transition. 



III. DEFECT PARTICLE 



In this section, we consider how condensation may arise when there is only a single 
particle of species B. The hop rates of the A particles are chosen to be 

u{n,0) = 1 and u{n,l)=p, (9) 

where p < 1, such that the A particles hop more slowly when the B particle is present 
(but they hop independently of n). Thus we view the B particle as a defect particle. The 
constraint on the hop rates (@)) then requires that the hop rate of the B particle is 

v{n, 0) = and v{n, 1) = . (10) 

Substituting these dynamics into 0, one finds that f{n,0) = 1 and f{n,l) = For 
this model, we can evaluate (0) as a finite sum (i.e. we can work directly in the canonical 
ensemble). Taking the B particle to be at site k, and summing over the L possibilities for 
k, yields 

L 

Zl,n,m = LY,S{J2^i-N)p-^'' (11) 

{m} 1=1 

This is the same normalisation (up to an overall factor of L) as that derived for the single 
species model with heterogeneous hop rates, in the case where particles hop from all sites 
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with rate 1 except for a single defect site from which particles hop with rate p \^ . Thus, using 

the results of one can identify two regimes: a low density phase, when pA < p/{l —p), 

and the system is in a fluid phase; and a high density phase, when pA > — p), and 

a Bose condensate forms at the site containing the defect particle. These regimes can be 

computed exactly by considering the normalisation (fT^ . and seeing that the sum may be 

dominated either by ~ C^fl)? in which case pa < p/i^ ~p), or by ~ 0(L), in which 

r 

case Pa > p/{^ ^ p) 8]. The condensate then contains a finite fraction of all the particles 
in the system — the remaining particles form a power law distributed background. It is a 
Bose condensate in the sense that the particles condense onto the site containing the defect 
particle (in the same way that in Bose condensation, the particles condense into the state 
of lowest energy: the equivalence is observed by identifying the site containing the defect 
particle in the zero-range process with the state of lowest energy in the Bose gas). 

A condensation transition of this kind persists as long as we have a finite number of 
(indistinguishable) defect particles: let's say there are M defect particles where in the limit 
L —>■ oo, we keep M fixed. Also, the hop rates for the A particles are u{n,m) = pm where 
m = 1, . . . , M. If the smallest of these hop rates is pi, then above a critical density of A 
particles, all the sites containing i particles of species B contain a finite fraction of all the 
particles of species A. But because there can only be a finite number of such sites, each of 
these sites must contain an infinite number of particles of species A — the condensate is 
distributed equally among the sites containing i particles of species B. 

The analysis of this section leads us to view the defect particle(s) as a disordered back- 
ground upon which A particles evolve. The special feature of the two species model is that 
this background may also evolve with prescribed dynamics. This is the case we consider in 
the following section, when the number of B particles is extensive. 

IV. FINITE DENSITIES OF BOTH SPECIES 

With the perspective of particle dynamics on an evolving disordered background, a case 
of particular interest in the two species zero-range process is when the evolution of one 
species, the B particles say, depends only on the number of particles of the other species at 
a site. Therefore we take v{n,m) = 1 + r(n) for m > 0, where r(n) is a general function of 
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n. Then from Q we deduce that f{n,m) is given by 

f{n, m) = [l + r(n)]-™s(n) , (13) 
where s{n) is another general function of n. We then use (j21) to infer the rates u{n,m): 

We assume in the following that r(n) is a monotonically decreasing function of n, and that 
in the limit n — oo, r{n) 0^. Inserting the form ()13|1 into (jZj), and performing the sum 
over m, yields 

F{z, y) = f2 s{n)z--^±f^ , (15) 

^^nz. y) = f2 ns{n)z--^±f^ , (16) 
oz ^-^ 1 + r n — V 

n=0 ' V / y 

= ; ';'">'„ . (17) 

ay ^ [1 + r(n) - y]2 

These equations determine z and given the densities pA and pb, via (jHl)- The radius of 
convergence of the sum over m is ?/ = 1 and we take, without loss of generality, the radius 
convergence of the sum over n to be 2; = 1. 

To analyse the possible transitions, we need to elucidate the behaviour of pA and pb when 
considered as a function of z and y. This will enable us to draw graphs of the dependences 
of Pa and pb ony, for fixed values of z, from which we can determine the densities for which 
the saddle point approximation remains valid. In particular, we wish to consider how pA 
and Pb change as z and y approach their radii of convergence — if pA or tends towards 
a finite value, then condensation ensues. To this end, we make the following observations: 

1. For fixed z, pA and pb are monotonically increasing functions of y. 

2. For 2 ^ 0, PA ^ 0. 

3. For y^O,pB^O. 

4. For z < 1, Pb is finite for all y (including y = 1). 



8 



We supplement these observations with the following three conditions on r{n) and s{n), 
which determine whether and pB converge to finite or infinite values when z and y 
approach their radii of convergence. For z —>■ 1 and ?/ < 1, if, as n — oo, 

> faster than 1/n , (18) 

1, if, as n ^ oo, 

^ faster than 1/n , (19) 
1, if, as n ^ oo, 

^0 faster than 1/n, (20) 

then Pb finite. These observations and conditions enumerate all possible ways that z 
and y approach their radii of convergence, and therefore all the possible circumstances in 
which condensation can occur in our model. The phase behaviour depends on which of the 
conditions (fTBj) to ()20|) are met and which are not. 

There are several possibilities, which we illustrate for a particular choice of r{n) and s{n), 
namely, for large n, 

s{n) ~ n~'^ , r(n) ~ cn~^ , (21) 

where h and c > are constants. Thus the asymptotic forms of the hop rates for large n are 
given by 

(1 - c/n^)"'(l + h/n) and v{n, m) ~ 1 + c/n . (22) 

We note that when c = the two species hop independently; in this case, the asymptotic 
.op ... oa,. . pa.«. .educe .o\ko. —ed i„ B to. single species zero-range 
process, where condensation was found above a critical density provided b > 2. 

With the choice (j2H), condition (|TH|l is satisfied if 6 > 2 and conditions (|T!Hl and are 
satisfied if 6 > 3. Therefore there are three cases to consider: 

Case 1: b < 2. 

In this case, none of the conditions (fTH|l to ((201) is met. The particular choice of rates 
studied in 10\ corresponds to this case. The dependences of pa and ps on y for fixed 
values of z are shown in Fig. Q where < Zi < Z2 < 1. Here, for a given p^, z 



ns{n) 

then Pa finite. For z 1 and y - 

ns{n) 
rin) 

then Pa —>■ finite. For z ^ 1 and y - 

s{n) 
r{nY 
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FIG. 1: Schematic dependences for case 1: b < 2, of the particle densities pA_ and pB for contours 
of fixed z and as a function of y. The dashed fine in the right hand graph illustrates how ps varies 
as a function of z and y given that pA is fixed (dashed line in left hand graph). 

must lie in the range zi < z < Z2. However, in this range, pb increases monotonically 
from pB = 0, where y = and z = Z2, to a. maximum value aX y = 1 and z = zi. If 
Pb exceeds this maximum then we can no longer solve the saddle point equations dHJ 
for both Pa and pB and the excess B particles condense onto a single site. Therefore 
whenever pb exceeds a pA-dependent maximum, the system is in a condensate phase. 
Otherwise the system is in a fluid phase. The critical line, given as a function of z for 
y = 1, is shown in Fig. [21 The explicit expression for the critical line is 

Pb = {1 + Pa)/c. (23) 

The condensate of B particles (which contains 0{L) particles) is induced by the 
distribution of A particles. In particular, at the site containing the B particle conden- 
sate, the A particles form a 'weak' condensate (which contains (9(L^/^) particles). To 
see this, note that the current of A particles must be finite, therefore u{n,m) must 
be finite at the condensate site. With the rates inferred from ()21|) . if m ^ oo then 
u{n, m) ^ unless we also have n ^ oo. Therefore taking n large one finds that 
u{n,m) ~ exp(— m/n^). Since this must be finite we must have m ~ at the con- 
densate site. Then, because m ~ 0{L), we must have n ~ 0{L^^'^). Away from the 
condensate site, the B particles form a power law distributed background and the A 
particles form an exponentially distributed background. 
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1/c PB 



FIG. 2: Phase diagram for case 1: b < 2. Phase I is a fluid phase; in phase II the B particles form 
a condensate sustained by a 'weak' condensate of A particles, as described in the text. cpA and cpB 
denote the numbers of particles contained in the condensates of A and B particles respectively. 

We can understand the weak condensate of A particles by considering a zero range 
process with a single defect site. Consider a single species of particles — the A particles 
— which hop with rate u{n) = 1 except at the defect site, where they hop with rate 
u{n) = exp(— (jfL/n^) {g is a constant). The hop rate from the defect site reflects the 
effect of the B particle condensate on the A particles in the two species model. At the 
defect site, u{n) — if n is small and so a condensate forms. But since the hop rate 
must remain finite, the condensate contains n ~ 0{L^^^) particles. Hence, from the 
perspective of the A particles, the condensate of B particles in the two species model 
plays the role of a defect site. 

Case 2: 2 < 6 < 3. 

In this case, only condition ()18|) is met: pA is finite for z = 1 provided y < 1. In this 
case, the dependences of pA and ps on y for fixed values of z are shown in Fig. 01 
again where < zi < Z2 < I. This time, imagine that the system is in the fluid phase, 
with densities pA and pb (and therefore values of z and y) corresponding to point A 
in Fig. ini Now, if we add more A particles to the system while keeping the density of 
B particles fixed, we find that we must increase z and decrease y as indicated by the 
dashed line. However, when z reaches 1, if we add more A particles to the system we 
can no longer solve the saddle point equations for z and y, therefore the A particles 
must undergo a transition from a fluid phase to a condensate phase. The critical line 
is given hj z = 1 — the critical density of A particles increases with increasing pB- We 
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FIG. 3: Schematic dependences for case 2: 2 < 6 < 3, of the particle densities pA and pB for 
contours of fixed z and as a function of y. See text for commentary. 




1/c PB 



FIG. 4: Phase diagram for case 2: 2 < 6 < 3. Phase I is a fluid phase; in phase II the B particles 
form a condensate sustained by a 'weak' condensate of A particles, as described in the text. A 
condensate of A particles and fluid of B particles form in phase III. (f)A and (j)B denote the numbers 
of particles contained in the condensates of A and B particles respectively. 

have not been able to find an explicit expression for this critical line. Also note that at 
?/ = 1, the system must undergo a transition between a fluid phase and a condensate 
of B particles as described in the previous case. Thus we deduce the phase diagram 
shown in Fig. EJ The two critical curves intersect when = oo and ps = oo. In phase 
III, the condensate of A particles exists on a power law distributed background of A 
particles while the B particles are exponentially distributed throughout the system. It 
is interesting to consider the sequence of transitions induced by increasing ps'- starting 
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PA 



cpA ~ 0{L) 



III 



IV 

^A ~ 0{L) 
4>B ~ OiL) 




II 
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FIG. 5: Phase diagrams for case 3: b > 3. In region I both species are in a fluid phase; in region II 
the B particles are in a condensate phase sustained by a weak condensate of A particles; in region 
III the A particles are in a condensate phase and the B particles form a fluid; in region IV both 
species are in a condensate phase. (pA and (pB denote the numbers of particles contained in the 
condensates of A and B particles respectively. 

from a point in phase III, the condensate of A particles is destroyed by increasing the 
density of B particles sufficiently, when the system enters the fluid phase I. Increasing 
Pb further leads the system to phase II where the B particles condense. 

Case 3: 6 > 3 

Here, all the conditions (fTH|). (fT^ and (^01) are satisfled. The arguments of the previous 
two cases apply but now, the critical curves given by 2; = 1 on the one hand and y = 1 
on the other intersect at flnite values of both pA and pB- Therefore when the A and 
B particle densities exceed their values given by 2 = 1 and y = I the system enters 
a phase where both species form a condensate at the same site. In this phase, the 
background distributions of A and B particles are both given by power laws. The 
phase diagram for this case is shown in Fig. El 

We have conflrmed the existence of the four phases presented in this section numerically. 
Exact expressions for P{n), the probability of finding exactly n A particles at a site, and 
P{m), the probability of finding exactly m B particles at a site, can be obtained in terms 
of the normalisation Zl^^^m- This normalisation satisfies an exact recursion equation j^l 



N 



M 




(24) 



71=0 m=0 
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FIG. 6: Log-log plot of on the left P{n) vs. n, and on the right P{m) vs. m, for systems of size 
L = 100 and 6 = 4 and c = 2. The circles correspond to densities pA = 1/2 = pb (fluid phase); 
the crosses correspond to densities pA = 1/2 and pB = 3 (condensate of B particles is sustained 
by a 'weak' condensate of A particles); the squares correspond to densities pA = ^ /'^ = Pb (both 
species condense); the diamonds correspond to densities /o^ = 3 and pB = 1/2 {A particles form a 
condensate). 

which is easily obtained from (0), and which can be iterated on a computer. Doing so, for 
systems up to size L = 100 with 6 = 4 and c = 2, yields the distributions shown in Fig. |^ 
Four phases are evident. The circles represent densities pa = 1/2 = ps and both species 
are in a fluid phase (phase I). The crosses represent densities pA = 1/2 and pB = 3, and 
the distributions are consistent with a condensate of B particles on a power law distributed 
background, and a 'weak' condensate of A particles containing n ~ C(Li/2) particles, on an 
exponentially distributed background (phase II). The diamonds represent densities pA = 3 
and pb = 1/2 and the A particles form a condensate on a power law distributed background, 
while the B particles form a fluid (phase III). The squares represent densities pa = and 
Pb = 7/2 and both species form condensates on power law distributed backgrounds (phase 
IV). 

It is possible to generalise the choice for r{n) in (j2H) to, for example, r(n) ~ cn~'^ where 
(i > is a constant. This does not lead to any phase diagrams topologically distinct from 
those already presented, although it does lead to the possibility of pb converging to a finite 
value as Pa oo, and vice versa, as z and y 1. Thus the phase diagrams in Figs. |21 
and 13 may be modified such that the phase boundaries tend toward a finite value of pb as 
Pa ^ oo, and vice versa also in the case of Fig. |3] Another feature of this generalised choice 
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for r(n) is that in the phase II, where the B particles condense, the accompanying 'weak' 
condensate of A particles contains a number of particles n ~ 0{L^^^^~^^^), as may be verified 
using the argument expressed in case 1. 



V. RELATION TO AHR MODEL 



In this section, we show that steady state of the two species zero-range process has a 
mapping on to the steady state of the AHR model. 

The AHR model, introduced in |[5;] , is a generalisation of the second-class particle system 
studied in It is defined on a ring of L + A + M sites, on which there are N + particles, 
M — particles, and L vacancies (which we represent by O's). The dynamics are defined by 
the processes 

+ 0^0+, with rate /3 , 
— — , with rate a , 
H — — h , with rate 1 , 

— h H — , with rate q , (25) 

where each exchange takes place between nearest neighbour sites. For a = j3 = 1 the model 
undergoes a transition between a disordered phase (g < 1) and a phase separated phase 
(g > 1) composed of a single domain of each species. The correspondence between the AHR 
model and the two species zero-range process may be observed in the following way. 

We define w{{Ti}) to be the steady state weight for the system to be in a configuration 
il}} = Ti, • • • , tl+n+m- The weights w{{Ti}) can be obtained using a matrix ansatz 



15|, that is, we write the particle configuration as a product of matrices 

{ti} = Xi-- ■ Xl^n+m , (26) 

where the matrix Xi is 

D iiT, = + , 

X = < E if Ti = - , 
A if Ti = . 

Then it can be shown that the steady state weights can be written in the form jl^ . Iisj ] 

«;({r,}) = Tr[Xi---XL+N+M], (27) 
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provided the matrices D, E and A satisfy the relations 

I3DA = A , (28) 
aAE = A , (29) 
DE - qED = D + E . (30) 

These relations are satisfied if we take A to be the projector where we employ a 

bra-ket notation to denote the left and right vectors {W\ and \V). With this notation, (f?fjl 
becomes 

w{{r^}) = {W\X,---Xl+n+m\V) , (31) 

and if the l-th vacancy is at site ki, then this can be written (choosing the normalisation 
(iy|V^) = 1 and using the invariance of the trace under cyclic permutations of the X's) 

L 

= n(^l^'^.+i ■ ■ ■ Xh^^-i\V) , (32) 

1=1 

i.e. the steady state weights assume a factorised form — one factor for each vacancy. Now, 
to make the connection to the two species zero range process, we define the matrix G^.m 
to be the sum over all permutations of products of n D's and m E^s. Also, we define 
P{{ni}; {mi}) to be the probability that, in between all pairs of vacancies / and / + 1, there 
are exactly rii + particles and rrii — particles. Hence 

L 

P{{ni}; {mi}) = Z^^^^,,l[{W\Gr,,m^V) , (33) 

1=1 

where Z^^^^m is a normalisation. is identical to (H)) if we make the identification 

f{n,m) = {W\Gn^m\V) , (34) 

and the normalisation Z^^^ m then is given by This establishes the mapping. 

Thus the steady state of the AHR model can be expressed in a form identical to the 
steady state of the two species zero range process if we identify the + particles with the A 
particles and the — particles with the B particles. The hop rates of the A and B particles, 
obtained by substituting (jMj) into are given by 

O^-'l^-'-l"* aud .(,,,n)=<^%^. (35) 



u[n,m] 



{W\Gr.,m\V) ' ' ' {W\Gn,m\V) 
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Note that because the mapping specifies f{n,m) (and not the hop rates) the hop rates 
are guaranteed to satisfy the constraint equation (jl)). Also, this is not a mapping for the 
microscopic dynamics — rather, it is a mapping between steady states which have the same 
form. 

The matrix elements (W^IGn^miy) are known exactly and assume different asymptotic 
forms depending on the values of a, (3, and q. Thus the phase behaviour of the AHR model 
is observed in the two species zero-range process by using these different forms to determine 
f{n,m) using (jS^)- We note that the resulting values of the hop rates are different to those 
studied earlier in this paper: the hop rates as determined via the mapping to the AHR 
model obey the symmetry u{n^m) = v{m^n) under the interchange a ^ (3. 



VI. CONCLUSION 



We have shown how the steady state of the two species zero-range process can undergo 
a number of condensation transitions. A single particle of one species was found to be able 
to induce condensation in the other above a critical density. Next, for finite densities of 
both species, we investigated a case where the hop rates of the two species were coupled in 
a nontrivial way. Three distinct condensate phases emerged and the conditions on the hop 
rates leading to such phases were presented for quite general rates. This generality suggests 
that the transition mechanisms are robust. 

There remain a number of outstanding questions. A more detailed understanding of the 
phase where the condensate is sustained by a 'weak' condensate of particles of the other 
species is desirable. It is also unclear whether there exist further couplings between the 
particle species which might lead to new transitions. This could require analysis of the 
model for dynamics which do not satisfy the constraint (Q; such investigation may also 
yield insight into the structure of the steady state when the factorised form does not hold. 
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